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ABSTRACT 
We propose a non-Gaussianity test for gravitational wave backgrounds by combining data streams of 
multiple detectors. This simple method allows us to check whether a detected background is "smooth" 
enough to be consistent with an inflation-type background, or is contaminated by individually unde- 
tectable weak burst signals. The proposed test would be quite useful for the Big Bang Observer or 
DECIGO whose primary target is a background from inflation at 0.1-lHz where gravitational wave 
bursts from supernovae of population III stars might become a troublesome foreground. 
Subject headings: gravitational waves — early universe — supernovae: general 



1. INTRODUCTION 

Since gravitational interaction is very weak, a gravitational wave (GW) background can serve as an invaluable fossil 

f-H from early universe with almost no scattering and absorption during its propagation (Maggiore 2000). By analyzing 

f-^ the background, we might obtain crucial information to understand physics at extremely high energy scale, e.g. the 

I inflation process. However, overlap of multiple signals in data streams of detectors would become a basic aspect of GW 

O ' astronomy, especially in the low frequency regime where number and variety of astrophysical sources would increase. 

^ I This is partly because GW detectors have omni-directional sensitivity. Therefore, for detecting a background from 

^ i early universe, it is essential to disentangle diversified signals contained in observed data. For example, we need to 

I i ' detect and subtract individual astrophysical signals such as chirping binaries, from data streams (see e.g. Arnaud et 

' al. 2007; Harms et al. 2008 for recent progress). 

y—i • The primary target for the Big Bang Observer (BBO, Phinney et al. 2003) and the DECihertz Interferometer 

J> ' Gravitational wave Observatory (DECIGO, Seto et al. 2001) is a background from inflation in the 0.1-lHz band where 

^~~l ', a relatively deep window of GWs has been expected to be opened. However, it was recently pointed out that burst-like 

l/^ ' GWs produced at supernova (SN) explosions of population-Ill (popIII) stars might become a problem for detecting 

^~~^ ' an inflation background around 0.1-lHz (Buonanno et al. 2004; Sandick et al. 2006; Suwa et al. 2007). While a 

" I , burst signal with a large amplitude might be handled in data analysis, a potential background composed by weak 

t — ■ undetectable bursts would be a formidable obstacle for adequately identifying an inflation background. Furthermore, 

f^ ■ in contrast to regular waveforms {e.g. from binaries), it is often difficult to accurately model waveforms for burst 

OO signals and/or to predict their characteristic amplitudes beforehand. To deal with such a situation, we certainly need 

^D ' to study GW backgrounds with efficient quantification methods beyond the traditional simple measure ^gw, the 

Kj^ energy density of backgrounds (see also Finn et al. 1999; Drasco & Flanagan 2002). In this paper we propose a 

• ^H , non-Gaussianity test for GW backgrounds in 0.1-lHz band, and discuss its prospects with BBO, setting GWs from 

j^ • the popIII SNe as our fiducial burst model. 

5-H This paper is organized as follows: in §2 we briefly describe data streams of BBO with summaries for notations. The 
basic idea behind our non-Gaussianity test is presented in §3. In §4 we make a statistical evaluation for our method. 

2. DATA STREAMS OF BBO 

For BBO, two LISA-type units (armlength L = 5 x lO^km) would be used to form a David-star-like configuration. 
From each unit we can obtain three data streams {A,E,T) and {A' ,E' ,T') (Prince et al. 2002; Corbin & Cornish 
2005; Seto 2006). Among these six data, T and T' modes are less sensitive to GWs, and we neglect them hereafter. 
The data streams A, E, A' and E' can be effectively regarded as responses of simple L-shaped detectors around their 
optimal frequency fopt ^ 0.3Hz with bandwidth Sf ^ J opt (Seto 2006). In this paper we only deal with quantities 
made from the following two pairs; A — A' or E — A' . The orientation of the arms of the former pair is aligned on the 
common detector plane of the two units, but the latter is misaligned by 45° on the plane. 

We model the data streams sx{t) {X — A,E and A') in terms of GW signal Hx and detector noises nx as 
sx(t) = Hx{t) + nx{t). For analyzing the noises nx, it is advantageous to work in the Fourier space. To clarify 
our main points, we discuss signal analysis in the optimal band with neglecting details of frequency dependence {e.g. 
replacing an integral J{- ■ ■)df with a product (• • O/opt ^ ^f)- I'^ practice, this situation is approximately realized by 
applying a band-pass filter. In order to take Fourier transformation, we decompose the data streams (total duration 
Toi,s) into short segments of a given period Tseg{^ fopt)t ^-^^ assign a label i\/(= 1, ■ • • ,Tobs/Tseg) for each segment 
with chronological order. Then we calculate 

sxM{f) = / sx{t)e^^'^'dt = Hxm{I) + nxM{f)- (1) 



The number of relevant Fourier modes in a segment is ^ TsegSf ^ Tsegfopt- In the Fourier transformations above, we 
imphcitly assumed to apply an adequate time window function to suppress leakage of underlying frequency components 
to nearby modes due to finiteness of Tseg ■ 

Hereafter, we only used the Fourier transformed quantities, assuming that the detector noises nx are stationary, 
Gaussian, and have identical spectrum SnU)- The assumption of Gaussian noises is not crucial for our method. For 
the relevant pairs {X,Y) = {A, A') or {E,A'), we also assume that their noises are independent (Phinney et al. 2003) 
with 

{nxM{f)nYL{f')*) ^ -^SMLSxYSff'TsegSNif), (2) 

where the notation (■ ■ •) represents to take an ensemble average. 

Next we discuss detectors' responses to incoming GWs. We expand the metric perturbation due to a GW background 
at a time t and position x by 

P=+,x / 

(c.c: complex conjugate) with three-dimensional wave vectors /, polarization bases ef^ and a definition / = |/|. 
Response of a detector X to an incident GW is characterized by the beam pattern function Fx{f, P)- In this paper, 
unless otherwise stated, we study simple L-shaped interferometers with the long-wave approximation. The explicit 
form of the function Fx is presented in the literature (Flanagan 1994; Allen & Romano 1997). For the background 
above, the signal HxAiif) in eq.([T|) is expressed as 

HxAlif)^ J2 E eM^T^Kf ■ XX + ftM)]hln{f,P)TsegFx{f,P) (4) 

with a shell-like three-dimensional frequency region Bf corresponding to the observed frequency /. Then we have 

{HxmUWymUY) = YSGWjifhxvTseg, (5) 

where the spectrum Saw if) for the background has dimension [Hz~^] as for the detector noise spectrum Sxif), and is 

written with the normahzed energy density ^gw by Scwif) = 327v"g f~^^GW (/) (-^o — 70km/sec/Mpc: the Hubble 
parameter). The overlap function jxrif) characterizes magnitude of common responses of two detectors {X,Y) to 
isotropic and unpolarized backgrounds (Flanagan 1994; Allen & Romano 1997). It is defined by 

jXYif) ^^ f dn{F+F+* + F^F^x*)^2../n.(a;„-a;,)^ ^g) 

o""" J Sphere 

We have ^eA' = due to geometrical symmetry and also ^aa' ~ 1 around the optimal band / ~ 0.3Hz of BBO (Seto 
2006). 

3. NON-GAUSSIANITY TEST 

If the signals HxKiif) in eq.(l4]) are made from superpositions of vast number of incoherent waves, they can be 
regarded as Gaussian variables from the central limit theorem. As a result, they are characterized only by second 
moments, and we have 

{{HxM{f)HYM{frf) = 2 {HxMimYMifY? cx ^Iy (7) 

from the properties such as {HxM{f)HxMif)) = 0. 

For a "coarse background" made by a relatively small number of freedom {e.g. popcorn noise due to super- 
nova bursts), the responses of detectors would be deviate from Gaussian. Therefore, using the BBO pair E-A' 
with the overlap function ^ea' = 0, we can check granularity of an isotropic background through the quantity 
(^[HEM{f)HA' M{f)*Y) that should vanish for a Gaussian background e.g. that generated at inflation. This is a 
key point in this paper. For detectors with a finite overlap 7xf ^ 0, we can generalize this method by introducing 
combinations such as {{HxM{f)HY M{f)*Y) —2 {HxM{f)HYMif)*) to subtract Gaussian component (similar to the 
definition of the Kurtosis parameter K4 for standard characterization of non-Gaussianity (see e.g. Racine & Cutler 
2007)). We do not pursue this direction further. But the underlying approach proposed in this paper would be appli- 
cable to a network of ground based detectors. For general detector configurations such as the LIGO- VIRGO pair, we 
need the subtraction scheme described above. 

For statistical analysis with BBO, we introduce the following two quantities made from the two pairs A — A' and 
E — A' respectively: 

C2=Y.Y. ^AM{f)sA'M{f)\ C^ = Y.Y. {^EM{f)sA'M{f)*f (8) 

M fdSf M f£Sf 

Here, the second summations X^fg^f ^'"^ fo'^ the Fourier modes (total number: TgegSf) within a segment, and the first 
ones X^A/ ^1'^ fo^ the segment label M — (1, • • • , Tots/Tseg). The combination C2 is used for traditional correlation 



analysis to measure flew, while C4 is our new probe for non-Gaussianity of a GW background. We evaluate their 
signal-to-noise ratios in the next section. 

In this paper we have set GWs from popIII SNe as our fiducial burst model. Here we comment on other models. 
Recently, several cosmological scenarios were proposed to produce intrinsically non-Gaussian GWs background in early 
universe. With a typical cosmological mechanism {e.g. through preheating phases), GW background is generated by 
causal processes, when wavelength were comparable to or smaller than horizon size in order of magnitude sense (see 
e.g. Easther & Lim 2006). Therefore, even if the generated GWs have a correlation structure, we have a large number 
(> {fHQ^y ^ 10'^^(//lHz)^) of independent emission regions for GWs currently observed at frequency /, and it 
would be difficult to directly probe the intrinsic non-Gaussianity for these typical models with BBC, due to the central 
limit theorem. But a background made by sparse cosmological events, such as GW bursts from cusps of cosmic strings 
might be an interesting target (Damour & Vilcnkin 2005). 

4. BURST BACKGROUND 
4.1. Derivation of Formulas 

In this section we analyze a GW background made by a superposition of burst events from single-species sources 
that have an event rate R and a characteristic duration time Td for GW signal in the optimal band. We start our 
discussion in a somewhat general way, and derive useful expressions for C2 and C4. Our fiducial model (GWs popIII 
SNe) will be examined later. 

We assume a smooth spectral profile A{f)^ {A > 0) for the GW emission throughout a burst event, and do not deal 
with regular waveforms with sharp frequency structures {e.g. monochromatic waves or chirping waves). For simplicity, 
the emitted wave pattern observed at the Earth is assumed to have an axisymmetric profile described by 

(/i+,/ix)cx(Q+(/),ax(/)) (9) 

(/: inclination angle, a+.x: complex number) in the coordinate system that is apparently symmetric to the geometry 
of the source. We also impose a normalization condition; i JJ^(|a+p -I- |ax P) sinldl — 1. 

For our probes C2 and C4, we need to evaluate the second- and forth-oder moments of the responses Hxm of 
two detectors induced by a burst event. Each response Hxm depends on the geometry of the source relative to the 
detector. The geometry is characterized by four angular parameters; the direction {9, 0) of the source on the sky and 
the orientation {1,1^) of its axis, and we have 

Hxm (X A{F+{0, 0, ^)a+(/) + F^ {0, 0, V)ax (/))• (10) 

Since these four geometric parameters are randomly distributed for extra-Galactic events, we define an averaging 
operator [• • ■]an with respect the direction and orientation of sources, and obtain 

[HxMHYM]an — -^^fXY ■ (H) 

5 

We define the ratio Q = [{Ham H"^, jy,j)'^]an/ {[Ham H'^, j^j]an)'^ for quantitative evaluation of the probe C4. Its numerator 
is explicitly given as 

[{HAMH\,M)Xn^^l (9{|a+|4-f |ax|^}-34|a+ax|' + 52Re[(a;ax)'])sin/d/. (12) 

Unless a weird cancellation occurs, the ratio Q becomes order of unity. At the end of this subsection, we will explicitly 
demonstrate this for our fiducial popIII SNe model. 

With a relation 7xx — 1 for a self correlation, eq. (fTTj) shows that the angular average of the response function is 
l/-\/5. Therefore, the characteristic signal-to-noise ratio SNRsst for individual burst with a single detector is evaluated 
as (see e.g. Segalis & Ori 2001; Sago et al. 2004) 

SNRBst ~ ^ ' , :,:' ■ 13 

51/25^/2 

Next we analyze the GW background formed by incoherent superposition of the bursts analyzed above, assuming 
that the typical signal-to- noise ratio SNRbsi is not larger than 0{\). As the expected number of events in a segment 
is RTseg, we get the background spectrum 

s -^^-^'^ (U) 

JGW,Bst — 7, U^J 

from expression (O and the corresponding mean magnitude (HxmHym) = RTseg x {A^) "ixY /^- Now we evaluate the 
signal-to-noise ratio SNRBst.C2 for the quantity C2 under the condition Tseg S> Td- If a segment M contains a burst 
event, its averaged contribution to C2 is (A^) 'fxyTsegSf/b. After incoherent superposition of TobsR{'> 1) events 
during observation time Tobs, we obtain the expectation value (signal strength) (C2) — (A^) IxvTsegSf /b x {TobsR)- 
Meanwhile, assuming independence of detector noises and the condition SN{fopt) > SGw{fopt), the rms fiuctuations 
for the product sams\ij^j are given by the detector noise spectrum as 2'^^''^2^^SNTseg- Here the additional factor 



2 ^ " is associated with the projection operation of data to the expected phase direction (usuaUy onto the real axis) in 
the complex plane (Seto 2006). Taking into account the total number of Fourier modes [Tobs/Tseg) x (TsegSf) — TobsSf, 
we obtain the rms fluctuations for C2 as 2'~^/'^SjsiTseg{Tobs5fY^'^, and its signal-to-noise ratio is given by 



^^i^(2T.,.5/)V2. (15) 



SNRBstfi2 = , „ 

5 ojv 

We can derive the same result for Tseg ^ Td- After replacing the product X(5/ with a frequency integral J df and using 
the spectrum R (.A^) = ^^^SGw,Bst^ the square value of this expression exactly matches with the standard formula for 
correlation analysis given in the literature (Flanagan 1994; Allen & Romano 1997). 
In the same manner, we can derive (C4) — (-4"*) QRTgegTobsSf /25 and obtain its signal-to-noise ratio as 

^NU 2y^A^)\Q\R{nbsSfr/' 

SNRBst,Ci = ^ — 7^^ — ^ (for Tseg > Td) (16) 

with an explicit dependence on T^eg in contrast to C2. Note that, here, independence of detector noises is an important 
requirement, but their Gaussianity is not essential. 
For Tseg ^ Td, we get the expectation value (C4) = (-4"*) QRT^^ TobsT^ <5//25 with its signal-to-noise ratio 

^NJ, 2^I^A^)\Q\R[T,bs5ff/^ 

SNRBst.Ci = ^ 25I\S2 (for Tseg < Td). (17) 

Eqs.dUl) and ^ are re-expressed as SNRBst.ci = 4Sfm!,J{T^,n,,\ ^^^Bst,c2-SNRl^^ with U = \Q\ (A^) / (A^)^ = 
0(1). With increasing the segment length Tseg from its minimum ~ fopt^ ^^ have a transition from (C4) oc T^^g to 
(C4) ex Tseg at the point Tseg ^ Td where the signal-to-noise ratio SNRBst,C4 also starts to decrease from a constant 
given in eq. (|17p due to dilution of power. Therefore, if the signal C4 is detectable at Tseg ^ fopt^ ^^ *^an estimate 
the duration time Td by identifying the transition. In the following, we set Tseg ~ fopt that will provide us with the 
maximum value of SNRBst.a for a burst background. 

Now we focus our discussion on a burst model with Td ^ f^Jf. {e.g. for our fiducial popIII SNe model at G.l-lHz). In 

this case, we have SNRBst,C4 — U-SNRBst,C2-SNR'^^^/A, and this relation is very suggestive. Even if individual burst 

1 /? 
has signal-to-noise ratio SNRbsi less than detection threshold, there is an amplification factor SNRBst,C2{oc T^^^ ) 

that increases with observational time Tabs and enable us to statistically study the bursts. 

With the parameters related to the bursts, we obtain SNRBst.a = {SGW,BstfU/R x MTT'^{2TobsSff'^5f/{2hS%). 
For a fixed background level SGw,Bst{o< R {A^)), the signal-to-noise ratio SNRBst,C4 decreases for a higher event rate 
R (corresponding to a smaller amplitude {A"^)). This is reasonable, considering that the background would become 
more Gaussian-like. If the bursts events are supposed to be the dominant sources of the total GW background around 

fopt and both (C2) oc R(A^) and (C4) oc R(^A?) U/Td are measured, we can roughly estimate the event rate R and 
the amplitude (-4^) separately, assuming U — 0{1). In addition to the estimated duration Td, these will be basic 
information to disclose the nature of the burst sources. 

As we commented earlier, the ratio Q = [{HEMH\,f^jY]an/i\HAMH\,i^,j\anY becomes order of unity for typical burst 
waveforms. Here we demonstrate this for our fiducial model: burst GWs from popIII SNe. In the band around 0.1-lHz, 
the emitted waves are dominated by memory effects caused by anisotropic neutrino emissions at supernova explosions 
(Buonanno et al. 2004), and we put the GW waveforms as (a+, a^ ) oc {sin^ /(I -I- |<7 — igsin^ /), 0} (Epstein 1978). 
The parameter q(< 0(1)) characterizes asphericity of the emission with q > for a polar enhancement and q < Q for 
an equatorial enhancement. Compared with the standard expression, we omitted an overall factor (oc q) that can be 
absorbed into the power A. For this model the ratio is given as Q = 0.99 [q = 1), Q = 1.02 {q — 0) and Q — 1.10 
[q = —1). We have limited our non-Gaussianity test only with the fixed pair E — A' . By combining results from 
geometrically different pairs, we can in principle, extract polarization information for burst signals. 

4.2. Prospects around O.TlHz 

In this subsection we specifically discuss prospects of our non-Gaussianity test for popIII SNe GW background with 
BBO. Since the anisotropics of neutrino emissions from popIII SNe or the formation rate of popIII stars are poorly 
understood, the amplitude of the popIII background is currently quite uncertain. With a parameter set for popIII SNe 
(redshift z ^ 15, emitted neutrino energy E^, ~ lO^^erg with mean isotropy {q) ^ 0.03 and event rate R ~ O.Olsec"^) 
extracted from Buonanno et al. (2004), we have ^gw ~ 4 x lO"""^^ at / '^ O.SHz. Meanwhile, for bursts characterized 
by the background level ^gw and their rate R, BBO has a sensitivity corresponding to 



for individual bursts, and 
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for the background with Tseg -Sf ^ 1, Q ^ 1 and U ^ \. Therefore, while identification of each burst might be difficult 
with small SNRbsi , our method has potential to discriminate whether a background once detected is smooth enough 
and consistent with inflation origin. 

An interesting question related to our non-Gaussianity test is whether we can separate smooth and burst contribu- 
tions for the total energy spectrum 0,gw- To estimate the latter component, we need the combination R (-4^)- With 
our approach based on a forth-order moment, we can obtain a different combination RQ {A^) — {Ri^A^))'^ x U/R. If 
the burst rate R is independently estimated e.g. with optical observation of popIII SNe, we can roughly estimate the 
burst component ex R(A^\ in the total spectrum ^gw by introducing a model parameter U — 0(1). 
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